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Abstract 



5— ( ' Exact diagonalizations with a minimally modified realistic force lead to detailed agree- 

ment with measured level schemes and electromagnetic transitions in ^^Ca, ^'^Sc, ^^Ti, ^^V, 
^^Cr and "^^Mn. Gamow- Teller strength functions are systematically calculated and repro- 
duce the data to within the standard quenching factor. Their fine structure indicates that 
fragmentation makes much strength unobservable. As a by-product, the calculations sug- 
gest a microscopic description of the onset of rotational motion. The spectroscopic quality 
of the results provides strong arguments in favour of the general validity of monopole cor- 
rected realistic forces, which is discussed. 



1 Introduction 



Exact diagonalizations in a full major oscillator shell are the privileged tools for spectroscopic 
studies up to A ^ 60. The total number of states — 2'^ with d = 12, 24 and 40 in the p, sd, 
and pf shells respectively — increases so fast that three generations of computers and computer 
codes have been necessary to move from n = 4 to n = 8 in the pf shell i.e. from four to eight 
valence particles, which is our subject. 

A peculiarity of the pf shell is that a minimally modified realistic interaction has been 
waiting for some 15 years to be tested in exact calculations with sufficiently large number of 
particles, and n = 8 happens to be the smallest for which success was practically guaranteed. 
In the test, spectra and electromagnetic transitions will be given due place but the emphasis 
will go to processes governed by spin operators: beta decays, (p, n) and (n, p) reactions. They 
are interesting — perhaps fascinating is a better word — on two counts. They demand a firm 
understanding of not simply a few, but very many levels of given J and they raise the problem 
of quenching of the Gamow-Teller strength. 

As by-products, the calculations provide clues on rotational motion and some helpful indi- 
cations about possible truncations of the spaces. The paper is arranged as follows. 

Section 2 contains the definition of the operators. Some preliminary comments on the 
interaction are made. 

In each of the following six sections, next to the name of the nucleus to which it is devoted, 
the title contains a comment directing attention to a point of interest. The one for ^^Ca is 
somewhat anomalous. 

In section 9, the evidence collected on GT strength is analyzed. Our calculations reproduce 
the data once the standard quenching factor, (0.77)^, is adopted. The fine structure of the 
strength function indicates that fragmentation makes impossible the observation of many peaks. 
Several experimental checks are suggested. Minor discrepancies with the data are attributed 
to small uncertainties in the a ■ a and err • ar contributions to the force. 

In section 10 we examine the following question: 
Why, in the sd shell, phenomenologically fitted matrix elements have been so far necessary to 
yield results of a quality comparable with the ones we obtain here with a minimally modified 
realistic interaction? The short answer is that monopole corrected realistic forces are valid in 
general, but the fact is easier to detect in the pf shell. 

Section 11 contains a brief note on binding energies. In section 12 we conclude. 

The rest of the introduction is devoted to a point of notation, a review of previous work 
(perhaps not exhaustive enough, for which we apologize) and a word on the diagonalizations. 
Notations. Throughout the paper / stands for f'i/2 (except of course when we speak of the 
pf shell) and r, generically, for any or all of the other subshells {pi/2 P3/2 fb/2)- Spaces of the 
type 



■no 



_|_ 



■no- 



(1) 



represent possible truncations: Uq is different from zero if more tlian 8 neutrons are present 
and wlien t = n — uq we liave tlie full space (p/)" for A = 40 + n. 

Bibliographical note. The characteristic that makes the pf shell unique in the periodic table 
is that at t = we already obtain a very reasonable model space, as demonstrated in the /" 
case (i.e. uq = 0) by Ginocchio and French [Q and Mc CuUen, Bayman and Zamick (MBZ in 
what follows). The uq ^ nuclei are technically more demanding but the t = approximation 
is again excellent (Horie and Ogawa 0). 

The first systematic study of the truncation hierarchy was undertaken by Pasquini and 
Zuker |Q, |^ who found that t = 1 has beneficial effects, t = 2 may be dangerous and even 
nonsensical, while t = 3 restored sense in the only non trivial case tractable at the time (^^Ni). 
This seems to be generally a fair approximation, and the topic will be discussed as we proceed. 

Much work on 7^ nuclei with t = 1 and 2 spaces was done by the groups in Utrecht 
(Glaudemans, Van Hees and Mooy P, 0, § and Tokyo (Horie, Muto, Yokohama ||^, [1^, 

For the Ca isotopes high t and even full space diagonalizations have been possible (Halbert, 
Wildenthal and Mc Grory 0, g). 

Exact calculations for n < 4 are due to Mc Grory [|14|, and some n = 5 nuclei were studied 
by Cole [|l^ and by Richter, Van der Merwe, Julies and Brown in a paper in which two sets of 
interaction matrix elements are constructed . 

For n = 6, 7 and 8 the only exact results reported so far are those of the authors on the 
magnetic properties of the Ti isotopes IT^, and the double /3 decay calculations of ''^Ca by 
Ogawa and Horie []T8[, the authors and Engel, Haxton and Vogel The MBZ model has 
been reviewed in an appendix to Pasquini's thesis ^ and by Kutschera, Brown and Ogawa ||21[ 



Its success suggested the implementation of a perturbative treatment in the full pf shell by 
Poves and Zuker [0, g. 

All the experimental results for which no explicit credit is given come from the recent com- 
pilation of Burrows for A = 48 . 



^^Ca 


^^Sc 


48Ti 


48Y 


48c^ 


m-scheme 12 022 


139 046 


634 744 


1489 168 


1963461 


J T 4 4 


4 3 


4 2 


5 1 


4 


1755 


17166 


63 757 


106 225 


58 219 



Table 1: m-scheme and maximal JT dimensions in the full pf shell 
The diagonalizations are performed in the m-scheme using a fast implementation of the 



Lanczos algorithm through the code ANTOINE |25]. Some details may be found in ref. 126]. The 



strength functions are obtained through Whitehead's prescription explained and illustrated 
in refs. [|l^, |19[ (and section 9). The m-scheme and maximal JT dimensions of the nuclei 
analyzed are given in table 1. To the best of our knowledge they are the largest attained so far. 
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2 The Interaction and other Operators 



The most interesting result of refs. [|, |] was that the spectroscopic catastrophes generated by 

in some nuclei, could be cured by the simple modification (KB' 



the Kuo-Brown interaction 
in I, I). 



19^(KB1) = V^,iKB) - (-)^300keV 
where V^^ are the centroids, defined for any two shells by 



Ej(2J + l) 



(2) 



(3) 



where the sums run over Pauli allowed values if i = j, and Wfjjj are two body matrix elements. 
For i,j = r, no defects can be detected until much higher in the pf region. On the contrary, 
the calculations are quite sensitive to changes in ^jjff but the only ones that are compulsory 
affect the centroids and it is the binding energies that are sensitive to them: 



V^f{KBl) = Vj/(KB) - 350 keV 
V}f(KBl) = y//KB) - 110 keV 



(4) 



The interaction we use in the paper, KB3, was defined in ref. |23| as 



Wf^ff(KB3) = W^^f f{KBl) - 300 keV for J = 1, 3 

W}jff{KB3) = 1^^21^/KBl) - 200keV. (5) 

while the other matrix elements are modified so as to keep the centroids (4). 

It should be understood that the minimal interaction is KBl in that the bad behaviour of 
the centroids refiects the bad saturation properties of the realistic potentials: if we do not accept 
corrections to the centroids we have no realistic interaction. Compared with the statements 
in eqs. (2) and (4), eq. (5) is very small talk that could just as well be ignored. However, to 
indulge in it is of some interest, as will become apparent in section 10. 

In what follows, and unless specified otherwise, we use 

• harmonic oscillator wave functions with b = 1.93 fm 

• bare electromagnetic factors in Ml transitions; effective charges of 1.5 e for protons and 
0.5 e for neutrons in the electric quadrupole transitions and moments. 

• Gamow- Teller (GT) strength defined through 

5(G^) = .^(ar)^ (^-) = ^^^=T^' (6) 
3 



where the matrix element is reduced with respect to the spin operator only (Racah con- 
vention |2^) and k is the axial to vector ratio for GT decays. 

= {9A/gv)eS = 0-77{gA/gv)hare = 0.963(7). (7) 



for the Fermi decays we have 

B{F) = (r)^ 

half-lives, t, are found through 



(fllE.^lli) 



(8) 



if A + nt 



6170 ±4 



(9) 



Uv/fA)B{F) + B{GT) 
We follow ref. |^ in the calculation of the and fv integrals and ref. for f^. The 
experimental energies are used. 



3 ^^Ca ERRATUM 

In fig. ^ we compare calculated and experimental levels. Except for the lowest 2+ state the 
agreement is good and the first excited 0"*" is certainly an intruder. 

The calculated Ml strength is found in a triplet (not shown) at excitation energies of 
9.82, 10.06 and 10.23 MeV with B{Ml) values of 0.39, 6.17 and 1.81 /x^, respectively, which 
nearly exhaust the sum rule of 8.96 /i^. The total observed strength between 7.5 and 12.5 
MeV is 5.2 ±0.5 yU^. It is dominated by the majestic peak at 10.2 MeV (3.9 and otherwise 
fragmented among some twenty states; below 11.7 MeV there are 14 peaks where the calculation 
only produce 8. The observed to calculated ratio 5.2/8.96 = (0.76)^, is very much the standard 
value for spin-like operators. 

The E2 rates 



B{E2,4+ 2+) = 2.65 e^fm^ 

B{E2,2+ 0+) = 10.2 e^fm^ (10) 
agree reasonably with the experimental values 

5(^2,4+ ^ 2+)exp = 1.89 e^fm^ 

B{E2X 0+)exp = 17.2 e^fm^ (11) 

but definitely suggest that something is missing in a strict Qhuo calculation. 

In ref. [|I9| we studied 2z/ double (3 decay of "^^Ca and calculated the strength functions for 
the associated processes '^^Ca(p, n)^^Sc and '^^Ti(n,p)'^^Sc, and for the latter we have to offer 
the following: 

ERRATUM. The total ^^Ti(n, p) strength in ref. misses a factor 3/2. Hence the ^^Ca 

1/2 

4 



2i/ double (3 decay half-life has to be multiplied by a factor 2/3 to yield T1/2 = 3.7 ■ 10^^ yr. 



Energy Levels of 48Ca 
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Figure 1: Experimental and theoretical energy levels of Ca 



48 



Sc The famous seven 



The J = 1 - 7, T = 3, multiplet of 7^" states in ^^Sc can be related to 7^" grates in ^^Sc 
through a Racah coefficient (the Pandya-Talmi transformation). The operation is successful 
enough to have become a textbook example. There are some discrepancies that are removed 
by a t=l truncation, hailed in ref. P as a triumph of the realistic forces. These levels change 

I and then to the exact results, but they are 



very little in going from t=l to perturbative |22 



quite sensitive to changes in the VFj/// Kicitrix elements. In all the other nuclei the situation 
is reversed and figure provides the only example of an exact calculation that does not bring 
an improvement over the approximate ones. Note however that the agreement with data is 
definitely good, and extends to three levels at around 2 MeV that do not belong to the multiplet. 
Below 2.5 MeV there are a couple of 2+ states with no calculated counterparts, i.e. intruders. A 
more complete view of the density of intruders comes from table |^, where we have listed the 9 
calculated l"*" levels below 6 MeV against 13 experimental candidates ||32|. Immediately above 
6 MeV, diagonalizations yield level spacings of 100 keV. The number of intruders will also grow 
fast and one to one identifications become meaningless because the number of levels that can 
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Energy Levels of 48 Sc 
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Figure 2: Experimental and theoretical energy levels of ^^Sc 
be observed becomes small fraction of those present. Note that even some of the calculated 



states in table 


2| may have escaped detection (e 


g. the one 


at 5.23 MeV). 








Exp 2.52 
Th 2.39 


2.98 3.06 3.16 3.26 3.71 
2.91 3.46 


4.00 4.1^ 
3.95 


i 4.32 4.68 
4.49 4.67 


4.78 
5.23 


5.45 
5.49 


5.74 
5.79 



Table 2: 1+ states in Sc. Energies in MeV. 

^^Sc decays to ^^Ti viathe doublet of 6+ states at 3.33 and 3.51 MeV. The measured half-life, 
log ft values and branching ratios are 



Ti/2 = 43.7 h 
log/t(6^) = 5.53 %/3 = 90.7% 
log /t(6^) = 6.01 %/3 = 9.3% 



(12) 



while the calculated ones read 



Ti/2 = 29.14 h 

log /t(6|) = 5.34 %/3 = 96% (13) 

log /t(6^) = 6.09 %f3 = 4% 

We have here a first example of the extreme sensitivity of the half-lives to effects that are 
bound to produce very minor changes in other properties that are satisfactorily described, such 
as those of the 6"*" doublet (see next section). 

5 ^^Ti Intrinsic states 

Experimentally this is the richest of the A = 48 nuclei. The lines in fig. ^ connect theoretical 
levels (to the right) with observed ones (to the left). The agreement, which would be perfect if 
all lines were horizontal, is nevertheless quite good: the rms deviation for the 33 excited states 
is of 110 keV. Dots correspond to the first intruders detected for a given J value. 



g 

8 

6 

>^ 
D) 

03 5 

LU 

o 

C3 3 
"o 

X 2 
LU 

1 




48Ti level scheme 



0+ 1+ 2+ 3+ 4+ 5+ 6+ 7+ 8+ 9+ 10+11+12+ 



Figure 3: Experimental and theoretical energy level of '^^Ti. The left side of the lines cor- 
responds to the experimental value, the right side corresponds to the theoretical one. The 
isolated points are intruders states not included in our valence space 
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Tables |^ and ^ contain the information on E2 and Ml transitions, to which we may add the 
magnetic moment of the first 2+ state, calculated to be /i(2+) = 0.43 fi^, lowish with respect 
of the measured fi{2+) = 0.86 (38) fiN and /i(2+) = 1.12 (22) H- 







EXP 


Up? 


/'(''Sc) 




0^ 


142 ± 6 


92.1 


52.9 


4^ 


2^ 


111 ± 20 


136.6 


65.8 






104 ± 104 


36.0 


45.2 






13 ± 2 


20.3 


6.6 




o| 


22 ± 6 


4.3 


0.2 


6^ 


4| 


52 ± 5 


54.8 


39.3 


6^ 


4f 


56 ± 25 


55.5 


0.7 


6^ 


4f 


75 ± 30 


34.1 


29.8 


8f 


6J 


<14.5 


20.2 


4.6 


8f 


6^ 


<51.8 


63.4 


41.4 


82+ 


6^ 


74 ± 40 


46.6 


32.4 


82+ 


6f 




10.7 


2.2 


10+ 


8^ 


<46.6 


46.7 


32.9 


9^ 


9^ 


<39.4 


11.8 


11.7 



Table 3: '^^Ti E2 transitions, B{E2) in e^ fm^. /®('^^Sc) is a t=0 calculation with Wffff matrix 
elements taken from the spectrum of ^^Sc. 

There are very few discrepancies between the exact calculations [fp)^, and the data and 
they are hardly significant. Remarkably enough this could also be said of the calculations, 
for the Ml transitions and to some extent for the B{E2) with the spectacular exception in 
B{E2, 6^2 4+). The exact results build enough quadrupole coherence with standard effective 
charges but the truly important difference with the t=0 calculations shows in table ^ where 
we have collected the intrinsic quadrupole moment Qq, extracted from the spectroscopic ones 
through 

_ (J+l)(2J + 3) 
- 3K^-J{J+1) 

setting K = 0. For the exact calculations, Qq is also extracted through the rotational model 
prescription 

B{E2,J J -2) = -^e^\{JK20\J -2,K)Ql. (15) 

167r 
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n \ / 




EXP 






ot 

^ L 


1+ 


0.50 ± 0.08 


0.54 


0.73 


Ot 

^ L 


z 


0.50 ± 0.08 


0.41 




Of 




0.80 ± 0.06 


0.42 




2^ 

Z 


i 


0.5 ± 0.1 


0.55 


1.63 




2t 


<0.01 


0.55 


1.63 


"1 


^t 


<0.06 


0.39 


0.82 




i 


1.0 ± 0.2 


1.78 


3.45 


2+ 


2| 

1 


0.05 ± 0.01 


0.39 


0.49 


Qt 


6^ 


6.7 ± 3.3 


2.32 


8.96 


6t 


6t 

^1 


0.9 ± 0.4 


0.06 


0.00 


'-'1 


6| 
"1 


>0.36 


0.41 


1.00 


5^ 


Qt 


>0.9 


0.27 


0.68 


7t 

' i 


^2 


0.09+^4 

— U.U4 


0.15 


0.30 


i 


6t 


0.14+^-?,^ 

— U.Ub 


0.33 


0.00 


8t 


8^ 
^1 


1.61 ± 0.70 


1.59 


3.38 


n 


8t 


U.0_o.2 


0.82 


0.00 


n 


&t 


-0.03 


0.12 


0.57 


9t 


8t 


>0.95 


0.82 


0.00 


9t 


8t 


>0.36 


0.57 


1.55 


8t 


7t 


0.31S:? 


0.54 


1.09 


8t 


St 


0.2tri 


0.25 


0.00 


10+ 


^t 


o.etli 


3.60 


7.08 


11^- 


10+ 


>0.36 


2.58 


5.37 


lit 


lot 


>0.36 


0.81 


0.00 


12t 


nt 


0.5 ± 0.2 


2.32 


3.46 



Table 4: ^^Ti Ml transitions, B{M1) in /x^. See also caption to table 3. 
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21 


50.4 


-12.25 


39.17 


68 


4^ 


32.18 


-11.83 


18.47 


69 


6^ 


-44.25 


13 


13.99 




6+ 


42 


4 


-25.18 


42 


81 


32.54 


1.67 


31.52 


44 


10^ 


46.69 


13.34 


30.29 


37 




35.78 


28.58 


31.64 






Qo(2 


J cxp 


62 ±3 





Table 5: ^^Ti intrinsic quadrupole moments of the Yrast states in e fm^. (/p)^ means Qq 
extracted from eq. (14). means Qo extracted from eq. (15). 



The wrong sign for the quadrupole moment is a characteristic of the t = calculations. The 
exact results check well with the known experimental value for the 2+ but they do much more: 
the curious mixture of signs and sizes coming out of becomes a fairly large and constant 
number for (fp)^. 

Although we cannot speak of rotational motion — which demands a truly constant Qo and a 
J(J+1) spectrum — we are certainly in the presence of a well defined prolate intrinsic structure. 

This build-up of quadrupole coherence is almost entirely due to mixing with the p^/2 orbit, 
as seen from the results for the (/7/2P3/2)^ space. 

The strength function for ^'^Ti(n,p)^^Sc can be found in ref. |]I9| (remember the erratum 
though. . . ) and ref. fl^ contains a study of the (p,p') and (e, e') processes and an analysis of 
the orbital, spin and Ml strengths. The missing piece of information, ^^Ti(p, n)^^V, is found 
in figs. ^ and ^, for t = 1 and exact calculations. They show similarity in gross structure: some 
low energy peaks, a resonance-like middle region and T = 2 satellite strength higher up. In 
details they differ mainly in the position of the resonance, shifted down by some 2 MeV in the 
exact case. The total strength for the full space = 13.263 k^, combined with = 1.263 
obtained for '^^Ti(n,p)^^Sc, satisfies the sum rule for the GT operator 



B{GT-) - B{GT+) =S- -S+ = ?,{N - Z) (16) 

The mirror (3^ decay of ^^Fe to ^^Mn has Qp = 11.1 MeV, which covers a large fraction 
of the strength in fig. |^. Experimentally, a comparison of the (3^ and {p, n) processes seems 
possible. It would certainly be welcome. 
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48Ti (p,n) 48V 




2 4 6 8 1012141618 20 22 

Energy (MeV) 

Figure 4: ^^Ti — > ^^V strength function, t = 1 calculation 



48Ti (p,n) 48V 




2 4 6 8 1012141618 20 22 
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Figure 5: ^^Ti — > *^^V strength function. Full fp-shell calculation 
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48 



V The interaction 



In fig. P we have plotted separately all levels with J <1 below 2.5 MeV and the high spin 
ones. To the repetitive claim that the agreement is excellent we may add a mitigating (i.e. 
possibly interesting) remark: it is next to impossible to observe the calculated 10"^ and 12^ 
states because of the screening provided by their close 11^ and IS"*" neighbours. 



Energy Levels of 48V 
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Figure 6: Experimental and theoretical energy levels of ^^V 

The Ml and E2 transitions in table ^ show again excellent agreement. For the ground 
state, /i(4+) = 1.90 /iAT against a measured /i(4+) = 2.01 (1) /iat, while /i(2+) = 0.51 /iat against 
two measures /i(2+) = 0.28 (10) /i^v and 0.444 (16) /i^ 

The Gamow- Teller /3+ decay to ^^Ti can reach J = 3, 4 and 5 daughter levels and three 
of them are fed as shown in table ^ The resulting half life T1/2 = 8.85 d is almost half the 
observed one, T1/2 = 15.97 d. The fraction of the total strength responsible for the decay is 
very small (0.4%). Therefore, the discrepancy in the half-lives could be cured by reducing by 
half the height of the lowest two bumps in figure ^, which is hardly a change in the overall 
picture. 
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59 


5+ 


42^ 


>176(124) 


41.0 


2t 


2f 


>1.3(19) 


10.7 






S(Ml)exp 


5(Ml)th 




2^ 


>0.027 


3.12 




4r 


0.081(14) 


0.188 


42^ 




0.045(9) 


0.032 


4^ 


4i^ 


0.0084(9) 


0.0079 


6^ 


5^ 


0.027(5) 


0.027 



Table 6: ^^V electromagnetic transitions, B{E2) in e^ fm^ and B{M1) in /x^ 



Jf 


B{GT),^ 


B{GT)exp 


4^ 


7.9 10-3 


4.1 10-3 


3^ 


0.5 10-3 


1.7 10-3 


42+ 


2.2 10-3 


4.1 10-3 



Table 7: (4+) ^ ^^Ti (Jf) beta decay 
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48V 48Ti 
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Figure 7: Strength function for the V beta decay. Full pf shell calculation 



Odd-odd nuclei provide a good test for the interactions and ^^V offers a good example 
of a general trend: perturbative calculations are quite good but the exact results come 
definitely closer to the data. This systematic improvement clearly indicates that the interaction 
is excellent. The one exception to the trend comes in ^^Sc, and is related to the choice of the 
KB3 (rather than KBl) interaction. As far as A = 48 nuclei are concerned this is an almost 
irrelevant problem, but there is a conceptual point that will be made clear in section 10. 
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7 ^^Cr Rotational motion 

Some levels in ^*Cr are populated in the decay of ^^Mn and will be discussed in the next section. 
Beyond that, little is known of the spectrum, except the Yrast line and before we come to it, 
we go through fig. |^ for the strength function. Only the first 1"*" state is seen through the Q13 
window. The calculated half-life is T1/2 = 21.8 h and the observed one T1/2 = 21.56 (3) h. 

48Cr^ 48V 
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CM 
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0.0 [ — , 



Figure 8: '^^^Cr ^*^V strength function. Full pf shell calculation 

Table ^ and fig. |^ collect the information about the Yrast levels. The only complaint we 
may have with the data is the B{E2, 8^ — 6^) value but the calculations are telling us much 
more than they agree with observations. 



J B{E2, J —> J — 2) Qspcc Qo 

exp. th. from Qspec from B{E2) 

2 321 ± 41 228 -29.5 103 107 

4 259 ± 83 312 -39.2 108 105 

6 >161 311 -39.7 99 100 

8 67 ± 23 285 -38.9 93 93 

10 >35 201 -22.5 52 77 

12 146 -5.3 12 



Table 8: Electromagnetic properties of the Yrast band of ^^Cr, B{E2) in e^ fm'', Q in e fm^ 
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Figure 9: Predicted and experimental Yrast band of Cr 

According to the criteria of constancy for Qo and J{J + 1) spacings defined in section 5 we 
cannot speak of a good rotor yet, but we are coming close to it. The J{J + 1) behaviour is at 
best incipient but the constancy of Qq is quite convincing. 

Two questions come naturally: what is the mechanism?, can we get better rotors? 

Table ^ suggests that it is the mixing of fj/2 and P3/2 orbits that is at the origin of a well 
defined prolate intrinsic state. If we follow this hint for any Aj = 2 sequences, it is easy 
to produce good rotors. The evidence is presented in table ^ It contains more information 
than necessary to make the point and should be viewed as introductory material for a future 
communication. Here we limit ourselves to a few comments. 

We use the bare KLS force |]3^ (all realistic interactions are very similar |^, ^] and 
for this one we have a code )• "^^^ oscillator frequency is hu = 9, the spacings between 



levels are of 1 MeV (higher j orbits coming below) and Qo{J) is defined in eq. (14). 

Under these conditions we obtain truly good rotors. Their stability has been tested under 
three types of variation. 



i) To affect all T = 1 matrix elements by a dilation factor A ~ 1.5 mocks surprisingly well 
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core polarization effects in a major sliell [35, Rotational behaviour persists and may 
be even emphasized. 

ii) For splittings between levels of 2MeV the rotational features are eroded but E4/E2 > 3 
in all cases. 

iii) Pairing is the most efficient enemy of good rotors. To go from KLS (bare) to KB (renor- 
malized), A = 1.5) (from i)) is sufficient, except for Wjj^j which should be doubled. Then 
E^^l E2 ~ 2.5, not far from what is seen in fig. ^ 
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0.992 
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0.990 


7 1 
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11.56 


0.955 
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0.995 
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0.979 


12.72 


0.979 


12 1 



Table 9: Excitation energies E (in MeV) and intrinsic quadrupole moments Qo (in fm^) for the 
Yrast states of several Aj = 2 configurations. The notation for the orbits is l2j. 

The relevance of these results to real rotors comes directly from Nilsson diagrams which tell 
us that the onset of deformation at ~ 60, 90 and 136 corresponds to the sudden promotion of 
two pairs of neutrons (u) and two pairs of protons (tt) to K = 1/2, 3/2 orbits originating in the 
spherical ones i^hii/2/n'gQ/2 z^^i3/2/7r/iii/2 and z^ji5/2/7rzi3/2 respectively. Naturally we identify 
the intrinsic state of the eight particles with the one coming out of the diagonalizations, which 
are possible right now for the first of the three regions, i.e. for the space (/^h/tPs)^ (ggd^si)^ 
(m-scheme dimension 2-10^). In table ^ we also show the results for {ggd^Si)^ (frPs)^, which 
we expect to become active for Z > 20, N ^ 40. The (ggd^Si)^ group is very likely responsible 
for rotational behaviour in the A = 80 region and may be even a candidate for "superbands" 
at lower masses. 
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8 ^^Mn Truncat ions and GT strength 

Spectroscopically, ^^Mn is identical to ^^V (to within Coulomb effects). Its decay to ^^Cr covers 



a non negligible fraction of the strength function |3g, |3^ . Since this process, and similar ones 
in the region, have been analyzed so far with t = 1 calculations, we are going to compare them 
with t = 3 and exact ones. A digression may be of use. 

In a decay, the parent is basically an /" state. Some daughters are also of /" type but most 
are /"~^r. In a t = 1 calculation, both configurations are present but is allowed to mix with 
jn-i^ through Wfffr matrix elements while f"~^r is not allowed to go to f^-'^r'^. At the t = 2 
level, pairing (i.e. Wffrr matrix elements) comes in and pushes /" down through mixing with 
while f^~^r cannot benefit from a similar push from f^-^r^. It is only at the t = 3 level 
that both /" and /"~^r states can be treated on equal footing. Now that t = 3 calculations 
are possible up to A = 58, it is useful to check their validity. 

The results for the half life of ^^Mn are collected in table We see that the t = 3 truncation 
is quite acceptable. For the total half-life the t = 1 number is not too bad if we remember 
that full calculations are not always as accurate as they are here. However this relatively good 
performance owes much to a strong Fermi branch: the partial GT value for t = 1 is simply 



bad. The last line of table |T0| gives an idea of the consequences of using calculated instead of 



observed energies for the transitions. 



half-life t = 1 t = 3 full exp 

Total 99 ms 133 ms 142 ms 158(22) ms 

Partial Fermi 244 ms 244 ms 244 ms 275(20) ms 

Partial GT 107 ms 292 ms 340 ms 372(20) ms 

Total (theoretical energies) 41 ms 80 ms 116 ms 



Table 10: ^^Mn half-lives with different truncations compared with the experimental results 



Next we move to the strength functions in figures [10], |TT| and 0, where the data show 
as dots. Globally the three calculations agree in that there is a resonance at ~ 13 MeV and 
structure at ~ 5 MeV. Beyond that, t = 3 provides a fairly accurate view of the exact shape 
while t = 1 does not. To give an idea of the energetics: the 4"*" T = 1 analog (IAS) of the ^^Mn 
ground state is at 5.79 MeV in ^^Cr, t = 1 puts it at 3.64 MeV, t = 3 at 4.77 MeV and the exact 
value is 5.38 MeV (not very good by our standards, see later in section 11). 

Now we come to the data and concentrate on the exact calculation. Figure |T3| contain very 



much the same information as figure [T^ but instead of (crr)^ we represent B{GT), affected by 



the (0.77)^ quenching factor for the calculated numbers. Furthermore, between 6 and 8.5 MeV 
we have eliminated among these the peaks that fall below the observation threshold (shown as 
a dashed curve) 
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Figure 10: ^^Mn — > ^^Cr strength function, t — 1 calculation 
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Figure 11: ^^Mn — > ^^Cr strength function, t — 3 calculation 
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Figure 12: ^^Mn — > ^^Cr strength function. Full calculation 
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Figure 13: Blow up of the low energy part of figure 12. Here we plot B{GT). Theoretical B{GT) 
values include standard (0.77)^ quenching. Dotted lines show the observational thresholds [39]. 
Bins of 250 keV for the calculations (vs. 500 keV in figs. 10-11). 

At low energy, the peak at 2 MeV and the cluster of states centered at 5 MeV are very well 
positioned in the calculations that yield a strength of 0.504 in the interval [0,5.75] MeV against 
an observed 0.485. Between 0.75 and 8.75 MeV nothing is seen experimentally, and the only 
calculated strength above the sensitivity limit is located in two peaks in the [7.25,7.50] MeV 
bin. At these energies the calculated levels are not yet eigenstates of the system but doorways 
whose strength (0.405) will be fragmented (see next section). 

Above 8.75 MeV observation resumes through delayed protons yielding B{GT) = 1.073 in 
the [8.75,11.75] MeV interval versus the calculated 2.07 mostly found in two bins at [10.25,10.50] 
and [11.25,11.50] MeV. At these energies the density of levels is high and fragmentation will 
become important. Furthermore, half of the calculated strength is in the last bin, dangerously 
close to the abrupt rise in the sensitivity threshold. Therefore, before we conclude that exper- 
iments demand anomalous quenching (i.e.: beyond the standard value) we shall analyse more 
closely what is being calculated and what is being measured. 
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9 Quenching, shifting and diluting GT strength 



From all we have said about GT transitions, a broad trend emerges: low lying levels are very 
well positioned and within minor discrepancies have the observed GT strength once standard 
quenching is applied. The examples of good energetics are particularly significant for the group 
around 5 MeV in ^^Cr and the l"*" levels in ^^Ti (fig. 0) and ^^Sc (table |]), that have experimental 
counterparts within 100 keV more often than not. The discrepancies are related to the shortish 
half-life of ^^Sc and ^^V, a slight lack of spin strength in the lowest states of ^^Ti (ref. [0) and 
— perhaps — with the tail of the resonant structure in ^^Cr discussed in the preceding section. 



48Mn^ 48Cr 
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Figure 14: Same as figure 12 but without any binning 

To cure the discrepancies we need a mechanism that may affect very slightly the overall GT 
distribution, without affecting the positions and other properties of the underlying levels, which 
are very well reproduced. Here, it is useful to remind how the strength function is obtained. 

First, we define a state |s) by acting on the parent, |s) = aT\i)S~^^'^ which exhausts by 
construction the sum rule S = {i\{aT)'^\i), and then evaluate the amplitude of |s) in each 



daughter state by doing Lanzos interactions using \s) as pivot [|17|. The first iteration produces 
the centroid and variance of the strength function. Eg = {s\H\s) and Vs = {s\H^ — Eg\s) 
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respectively. As the number of iterations v increases, the strength — originally concentrated at 
Es — is fragmented into v peaks. The lowest converge to exact eigenstates (with their exact 



share of the total sum rule) at the rate of roughly one every 6-10 iterations. Figure 113 is the 



high resolution version of fig. |T2| and shows the situation after 60 iterations for each of the 
JT values (J = 3,4,5; T = 0,1,2). It is only below approximately 6 MeV that we have a 
complete picture of the spectrum. Above, many thousand states are waiting to come and erode 
the strong peaks. 

Let us examine first the global properties associated with S and E^. and then turn to the 
consequences of local fragmentation. 

Quenching. Calculations always produce too much strength, that has to be reduced by a 



quenching factor, which is the stronger (i.e. smaller) the most drastic the truncation ||17|| . 
For exact calculations in one major shell (Ohuj), S has to be reduced by (0.77)^, which is 



very much the value demanded by the "violation" of the model independent sum rule ([16D 
There is very little we can do within a Ohu calculation to change this state of affairs. 



Shifting. Contrary to S, which depends on geometry (|TB]) and on overall properties of H, 
Eg may be significantly affected by small changes in the a ■ a and ar ■ err contributions 
to H. In ref. ^6 it is shown that these spin-spin terms are very strong, especially the 



second, and may differ from force to force by some 20% (the corresponding constants ei+o 



and Ci+i in ref. p6[ are related, but not trivially, to the Migdal parameters g and g'). 
Therefore, the mechanism to cure the small discrepancies we have mentioned may well 
come from modifications in these numbers, that would produce small overall shifts of the 
distribution, and nothing else. 

Now we return to the quenching problem. In view of discrepancies between [p, n) and 
data for ^"^Ca, the extraction of 5* from the former has been recently criticized 
The problem was compounded by the fact that calculations with Wildenthal's W interaction 
suggested that standard quenching did not seem necessary. This illusion was dispelled by 
Brown's analysis [^, showing that the interaction was at fault. It is very interesting to note 
that what is called 12. 5p in |^2|] is none other that KB, while CW — which gives the best 
results — is very basically a minimally modified KB that can be safely assimilated to KBl or 
KB3. 

Diluting. Brown's analysis contains another important hint: once normalized by the (0.77)^ 
factor, the CW calculations follow smoothly the data within the fi window, but then 
produce too much strength when compared with the (p, n) reaction. The suggestion from 
fig. is that in regions where the density of levels becomes high, fragmentation may be- 
come so strong that much strength will be so diluted as to be rendered unobservable. The 
difference between the A = 37 and A = 48 spectra is that in the former the effect is almost 
exclusively due to intruders, while in the latter the density of pf states is high enough 
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to produce substantial dilution by itself. Which brings us back to the problem of the 
amount of strength calculated between 8.75 and 11.75 MeV in ^'^Cr, double the measured 
value. Some shifting may be warranted to reduce the discrepancy, but fig. |T^ suggest very 
strongly that dilution be made responsible for it (i.e. for anomalous quenching) 

From all this, it follows that simultaneous measurements of (3^ and (p, n) strength are very 
much welcome in pairs of conjugate nuclei where the release energy is large, the Qhuj spectrum 
is dense and high quality calculations are feasible. We propose the following candidates: 

/3+ (p, n) 

46Qj, 46y _^ 46y 

47Mn ^ 47 47rji ^ 47V 

48Fe^48Mn 48^^ ^ 48v. 

To conclude: a theoretical understanding of standard quenching demands that we look 
at the full wavefunction and not only at its Qhuj components. Experimentally, what has to be 
explained is the disappearance of strength, i.e. standard and anomalous quenching (as observed 
in ^^Cr). Dilution will no doubt play a role in both, but the latter may be observed already by 
comparisons with Qhuo calculations. 



10 The validity of monopole modified realistic interac- 
tions 

To answer with some care the question raised in the introduction we review briefly the work 
related to monopole corrections. 

The flrst attempt to transpose the results of refs. |^ to the sd shell met with the problem 
that the interaction had to evolve from ^^O to ^°Ca. A linear evolution was assumed, but it 
was shown that the centroids followed more complicated laws demanding an excessive number 
of parameters . 



The solution came in ref. |3^, by adopting a hierarchy of centroids and noting that the 



realistic matrix elements depend on the harmonic oscillator frequency very much as 



w^/lH = -w^/lM, (17) 

thus displacing the problem of evolution of H to one of evolution of u. The classical estimate 
hu! = 40A~^/^ relies on flUing oscillator orbits and on adopting the r = tqA^^^ law for radii. 



which nuclei in the p and sd shells do not follow. Therefore it was decided to treat as a 
free parameter for each mass number and then check that the corresponding oscillator orbits 
reproduce the observed radii. This turned out to work very well and to produce very good 
spectroscopy in all regions where exact calculations could be done. 
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The method rehes on the rigorous decomposition of the full Hamiltonian as 7i = Tim + 
I-Lm-i where the monopole part Tim is responsible for saturation properties, while TIm contains 
all the other multipoles. Upon reduction to a model H,Hm is represented by Hm, which 
contains the binding energies of the closed shells, the single particle energies and the centroids. 
Everything else goes to Hm, which nevertheless depends on Tim through the orbits, in principle 
selfconsistently extracted from Tim [P^, 0- The program of minimal modifications now 



amounts to discard from the nucleon-nucleon potential the Tim part and accept all the rest, 
unless some irrefutable arguments show up for modifying something else. On the contrary Tim 
is assumed to be purely phenomenological and the information necessary to construct it comes 



mostly from masses and single particle energies |45 



The proof of the validity of the realistic TYa/ through shell model calculations depends on 
the quality of the monopole corrections. In regions of agitated radial behaviour we have to 
go beyond the oscillator approximation. In particular, the sd shell radii can be reproduced 
practically within error bars by Hartree-Fock calculations with Skyrme forces with orbital 
fillings extracted from the shell model wave functions ||4^. Obviously the ^5/2,^1/2 and ^3/2 
orbits are poorly reproduced by a single u, and obviously this makes a difference in the two 
body matrix elements (work is under progress on this problem). 

Therefore, in the sd shell — to match or better the energy agreements obtained with Wilden- 
thal's famous W (or USD) interaction ^ — we have to push a bit further the work of 



ref. 35 



When we move to the pf shell, no such efforts are, or were, necessary. Because of /7/2 
dominance, it is much easier to determine the centroids: the values are no issue, the crucial 
Vjj. ones can be read (almost) directly from single particle properties on ^^Ca and ^''Ni, and 
we are left with Vj^; only two numbers. Once we have good enough approximations for the 
centroids we can do shell model calculations to see how the rest of H (i.e. Hm) behaves. In 



ref. |]23| it was found that Hm behaves quite well, but much better in the second half of the 
region: A = 48 happens to be the border beyond which quite well becomes very well. Hence 
the remark in the second paragraph of the introduction. 

The trouble at the beginning of the region was attributed at the time to intruders, but 
now we know that radial behaviour must be granted its share. It is here that KB3 comes in. 
Although eq. ^ has only cosmetic effects, its origin is not cosmetic: it was meant to simulate 
necessary corrective action not related to the presence of intruders. It reflects the fact that in the 
neighborhood of ^°Ca the interaction works better if the spectrum of "^^Sc is better reproduced, 
while at A = 48 and after, the KB interaction — not very good in ^^Sc — needs no help, except 
in the centroids, as we have shown. 

Most probably this says something about radial behaviour at the beginning of the region. 



It is poorly known except in the Ca isotopes, where it is highly non trivial [49]. The indications 
are that when shell model calculations demand individual variations of matrix elements in 
going from nucleus to nucleus the most likely culprits are the single particle orbits. And the 
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indications are that at A = 48 the need for such considerations disappears. 

All this is to say that A = 48 happens to be a good place for a test of the interaction and 
KBl passes it with honours. It also amounts to say that it is not an accident but part of the 
general statement that monopole corrected realistic interactions are the ones that should be 
used in calculations. 

It is important to stress that monopole corrections are both necessary and sufficient. 

Doing less amounts to run great risks as illustrated by the T = 2 states in A = 56. Using 
KB leads to a very good looking spectrum in a t = 2 calculation of ^^Fe but ignores the 
fact that the force would produce nonsense for T = 0. (For T = 2, KBl or KB3 gives results 
that are even better looking than those of KB, and somewhat different). 

Doing more, under the form of extended fits of all matrix elements, has become unnecessary 
and may be misleading as illustrated by the problems encountered by W, the most famous 
of the fitted interactions (section 9, Still, a full explanation of its success in terms of 

monopole behaviour remains a challenge. 

11 Note on binding energies 

Our interaction overbinds all the A = 48 nuclei by about the same amount, indicating the 
need of small corrections in the centroids. The addition of an overall Vn{n — l)/2 term with 
V = 28 keV, leads to the binding energies in table [ll|, which have been Coulomb corrected by 
subtracting (n = vr + z/, vr = protons, z/ = neutrons) ^ 



Hcoui = K^vr(7r - l)/2 + V^^ttu + 7.2797r MeV 
Ktt = 0.300(50) MeV, 14^ = -0.065(15) MeV. (18) 





Ca 


Sc 


Ti 


V 


Cr 


Mn 


Fe 


Exp 
Th 


7.04 
7.10 


13.35 
13.35 


23.82 
23.79 


26.70 
26.80 


32.38 
32.17 


26.86 
26.80 


23.79 



Table 11: Nuclear two-body energies (MeV), A = 48 

For ^^Fe we obtain BE = 385.37 MeV, leading to a decay energy of 10.94 MeV for 
^^Fe — >■ ^^Mn. The most recent estimate of Wapstra, Audi and Hoekstra [ET| gives 



i?£'(^^Fe)=385.21 MeV. The precision of subtractions [^] can be checked through the posi- 
tion of the analog of ^^V in ^^Cr, experimentally at 5.79 MeV against the Coulomb corrected 



value of 5.68 MeV in table |TT[ As noted in section 8, for this state there is a 400 keV dis- 
crepancy between experiment and calculation, high by our standards. It probably indicates 
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the need of some more sophistication in the readjustment of monopole terms than the 28 keV 
we have suggested. Note that the theory-experiment differences in table |ll] are typical of the 
results throughout the paper. 



12 Conclusions 

Since much of the sections 9 and 10 was devoted to drawing conclusions about the two main 
problems addressed, we shall only sum them up: 

i) The detailed agreement with the data lends strong support to the claim that minimally 
monopole modified realistic forces are the natural and correct choice in structure calcu- 
lations. 

ii) The description of Gamow- Teller strength is quite consistent with the data, to within 
the standard quenching factor. The calculations strongly suggest that dilution due to 
fragmentation makes much of the strength unobservable. 

Two by-products emerge from the calculations. The first is mainly technical: 

iii) Truncations at the t = 3 level are reasonable in the lower part of the pf shell. In general, 
however, truncations are a dangerous tool and it would be preferable to replace them by 
some other approximation method. 

The second by-product is more interesting: 

iv) The calculations provide a very clean microscopic view of the notion of intrinsic states 
and of the conditions under which rotational motion sets in. 
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